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Orbitals from local RDMFT: Are they Kohn-Sham or Natural Orbitals? 
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Recently, an approximate theoretical framework was introduced, called local reduced density matrix 
functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1- 
RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron 
Schrodinger equation with a local potential. In the present work, we focus on the character of these 
optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained 
with the full minimization (without the extra condition) by contrasting them against the exact NOs 
and orbitals from a density functional calculation using the local density approximation (LDA). We 
hnd that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the 
full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of 
the description of strong static correlation, this finding opens the way to a mixed density/density 
matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the 
occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation 
at a cost only minimally higher than a density functional calculation. 


I. INTRODUCTION 


Reduced-density-matrix-functional theory (RDMFT) 
[I] is an alternative formulation of the many-electron 
problem where every ground-state property, inclnding 
the ground-state energy, is a functional of the one-body 
reduced density matrix (I-RDM). A main advantage 
compared to density functional theory (DFT) is that 
the electronic kinetic energy can be written explicitly in 
terms of the 1-RDM. Like in DFT, approximations in 
RDMFT can be cast in a form where all terms are sim¬ 
ple explicit functionals of the 1-RDM except a remaining 
unknown part of the electron-electron interaction term 
which can be also called exchange and correlation en¬ 
ergy. Contrary to DFT the exchange-correlation (xc) en¬ 
ergy does not contain kinetic-energy contributions since 
that part of the energy is treated exactly. Constructing 
functionals amounts to introducing approximate forms 
for the xc energy term in terms of the 1-RDM. Most ap¬ 
proximations in RDMFT are explicit functionals of - and 
are minimized in terms of - the natural orbitals (NOs), 
(/)j(r), and their occupation numbers, rij. Various differ¬ 
ent approximations for the total energy functional have 
become available over the past decades [2-17] which have 
proven to describe correctly such diverse properties as 
molecular dissociation [4-8] or band gaps[14, 18-20]. A 
major drawback of RDMFT, however, is the increase in 
computational cost compared to a DFT calculation which 
is mostly due to the optimization of the NOs. Contrary 


to DFT, the functional variation with respect to the or¬ 
bitals does not reduce to an iterative eigenvalue problem. 
A few techniques have been introduced to define effective 
Hamiltonian schemes with non-local potentials to obtain 
the natural orbitals [21-23]. Although some of these tech¬ 
niques reduce the computational cost substantially com¬ 
pared with the full minimization, orbital determination 
in RDMFT still remains a bottleneck when compared 
with DFT or even Hartree Fock methods. 

In an attempt to incorporate the merits of RDMFT 
functionals, like static correlation effects, in Kohn-Sham- 
like equations an alternative approach, called local- 
RDMFT, was introduced recently [24]. The main idea 
in this approach is to optimize the orbitals under the ad¬ 
ditional constraint that they are the eigenfunctions of a 
Hamiltonian containing only a kinetic term and a local 
scalar potential. Due to this strong constraint, local- 
RDMFT does not solve the full problem of orbital op¬ 
timization in RDMFT. Instead, the orbital optimization 
can be considered as being in the realm of DFT, and 
the orbitals satisfy single particle equations. The cost 
for one optimization of the orbitals in local-RDMFT is 
similar to the cost of a density functional calculation 
using an orbital functional via the optimized effective 
potential method. Local-RDMFT can potentially de¬ 
scribe static correlations, where most DFT approaches 
fail. More specifically, it was found that a qualitatively 
correct dissociation of H 2 and N 2 molecules is obtained. 
In particular, for H 2 , if the corresponding functional of 
the IRDM gives the correct physical picture of dissocia- 
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FIG. 1. The dependence of the occupation numbers of the Icr 
and Icr* spin-orbitals of H 2 as a function of the internuclear 
distance. Only functionals which yield a qualitatively correct 
dissociation are included (see text for details). 


tion in standard RDMFT, then this desirable feature is 
also preserved by local-RDMFT. 

In addition, the calculated energy eigenvalues of the 
Hamiltonian in local-RDMFT were shown to provide a 
useful energy spectrum for molecular systems. It was 
found that those energy eigenvalues are very good ap¬ 
proximations for the ionization potentials (IPs) of a set of 
atomic and molecular systems [24, 25] although a mathe¬ 
matical proof of their association with the IPs is lacking. 
In the present work, we refer to RDMFT with the full 
minimization of the NOs as full-RDMFT and use the 
prefix full also for specific functionals, e.g. full-Miiller, 
in order to distinguish from local-RDMFT for which the 
prefix local will be used. 

A given set of NOs cannot be obtained as the eigen¬ 
functions of a Hamiltonian with a local potential, as we 
know from their asymptotics. Hence, enforcing the addi¬ 
tional constraint is an additional approximation and the 
optimized local-RDMFT orbitals (LROs) are expected 
to be different from the natural orbitals. As a result, 
the question arises how much the LROs differ from the 
orbitals obtained from a full minimization, i.e. a min¬ 
imization of an approximate energy functional without 
the additional constraint, and from the exact NOs ob¬ 
tained with a very accurate multiconfigurational method. 
Since the constrained minimization scheme strongly re¬ 
sembles a density functional calculation one might expect 
the resulting LROs to be closer to the KS orbitals of a 
DFT calculation. This is certainly true for the asymp¬ 
totic behavior of the orbitals in finite systems. The exact 
NOs decay exponentially and if their occupation num¬ 
ber is fractional, which is generally the case, they will 
all share the same exponent for the decay, determined 
by the chemical potential of the system [26]. The KS 
orbitals for a finite system also decay exponentially but 
the exponent is determined by their KS eigenvalues which 
implies that the orbitals decay with different exponents 
unless they are energetically degenerate. The same holds 
true for the LROs, i.e. they decay with an exponent that 


FIG. 2. Helium 2s orbital along an axis through the position 
of the nuclei obtained with full- and local-RDMFT compared 
with the exact NO, and those obtained with LDA and RHF. 
The orbital obtained with either local-BBC3, local-Power or 
local-ML is very similar to the one given by local-Miiller. 


is determined by the associated energy eigenvalue. This 
has little effect on the total energy since the asymptotic 
region contains hardly any density and, hence, does not 
contribute significantly to the total energy. 

In the present work, we focus on the nature of LROs. 
More specifically, we compare the LROs from an approxi¬ 
mate RDMFT functional in the local framework with (a) 
the orbitals from a full minimization using the same ap¬ 
proximate functional, (b) very accurate approximations 
of the exact NOs obtained with multi-configuration self- 
consistent field (MCSCF) calculations, and (c) the KS 
orbitals from a DFT calculation using the local density 
approximation (LDA). This allows us to investigate if 
the orbitals from local-RDMFT are indeed closer to KS 
orbitals than to the exact NOs. We shall find that the 
LROs resemble the KS orbitals. Hence, our work pro¬ 
vides the motivation for further exploration of the possi¬ 
bility of hybrid DFT/RDMFT approaches, for example 
to obtain the orbitals from a flavor of DFT and the oc¬ 
cupation numbers through a minimization of an appro¬ 
priate RDMFT functional. 

The paper is structured as follows: in section H we re¬ 
view the basic ideas of local-RDMFT and the Schrddinger 
equation that the LROs have to satisfy. In section HI we 
then compare the shapes of the NOs from a full mini¬ 
mization with those from local-RDMFT, using different 
common approximations for the exchange-correlation en¬ 
ergy in RDMFT [2, 5, 6, 9, 27], for the helium atom and 
for the hydrogen molecule at different internuclear dis¬ 
tances. To complete the comparison we also include the 
almost exact NOs from a MCSCF calculation, the KS 
orbitals from a DFT calculation using the LDA, and the 
Hartree-Fock (HF) orbitals. We conclude the paper in 
section IV. 
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FIG. 3. H 2 Icr (left) and H 2 la* (right) orbitals at equilibrium plotted along the nuclear axis. For the Icr, all full- and local- 
RDMFT functionals employed yield the same orbital in the scale of the plot, therefore, only the local-Mliller orbital is shown 
which also almost coincides with the exact NO and that from RHF. For the Icr*, all local-RDMFT functionals employed yield 
almost the same orbital, while the full-RDMFT orbitals look slightly different from each other but are all close to the exact 
natural orbital. Only the full- and local-Mliller orbital is shown 


II. LOCAL RDMFT 


mentation of two additional constraints [28]: 


The theory of local-RDMFT was presented recently in 
[24], here we give a brief summary of the method and 
state the equations which are solved numerically. Com¬ 
monly, RDMFT functionals depend separately on the 
occupation numbers nj and the NOs, V^j(r). The to¬ 
tal energy is then minimized in a two-step process: first 
one minimizes with respect to the occupation numbers 
for a fixed set of NOs and then keeps those occupation 
numbers fixed while minimizing with respect to the NOs. 
These two steps are then iterated until convergence is 
achieved. While the minimization wrt. the occupation 
numbers is numerically inexpensive, the optimization of 
the NOs is quite costly. The local-RDMFT approach 
reduces the cost for this minimization significantly while 
still providing fractional occupation numbers through the 
usual optimization under Coleman’s iV-representability 
conditions. 

The central idea in local-RDMFT [24] is that the 
search for the set of LROs is restricted to the domain 
of orbital sets that satisfy single-particle equations with 
a local potential 


J dr prepir) = N - 1, (3) 

Prep(r) > 0 . (4) 

The two constraints are sufficient to correct the asymp¬ 
totic behavior of the effective local potential, thus healing 
to a large extend the self interaction problem expected to 
be present in the potential for many approximate func¬ 
tionals in local-RDMFT. 

The optimal ERD and the effective local potential can 
be obtained, similarly to the OEP method, by solving 
the integral equation [24]: 

J d^r' x(r, r') Prep(r') = b{r), (5) 


with 



x(x,y) 

|x-rlly-r'l ’ 


( 6 ) 


6(r) = 



&(x) 

[x-r] ■ 


(7) 


—— 1 14xt(r) + Frep(r) 


= ej(l)j{r). ( 1 ) 


The search for the electron-electron repulsive part Vrep(i") 
of the effective local potential (the analogue of the 
Hartree-exchange and correlation potential in the KS 
equations) is replaced by a search for the effective repul¬ 
sive density (ERD) Prep(r) whose electrostatic potential 
is Erep(r), i-e. 

V^Erep(l’) = -47rprep(r) . (2) 

Minimizing with respect to prep(r) allows for the imple- 


The response function x(r, r') and 6(r) are given by 

x(di-')= (8) 

pU) _ 

b{r)= -axc|^^)^*(r)0,(r), (9) 

with defined by 

= j d^r' F^l{r, r') . (10) 
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^'Hxc is the approximation for the electron-electron in¬ 
teraction energy, (j)j are the LROs and rij, and tj their 
corresponding occupation numbers and orbital energies 
(eigenvalues of the effective Hamiltonian). The two con¬ 
straints are incorporated with a Lagrange multiplier (3) 
and a penalty term (4) that introduces an energy cost for 
every point r where prep(r) becomes negative. 

Within local-RDMFT, the numerical cost for one min¬ 
imization wrt. the NOs reduces to the cost for a DFT 
calculation using an orbital functional. Apart from re¬ 
ducing the numerical cost, the local-RDMFT framework 
provides an energy eigenvalue spectrum connected to the 
LROs. We found that the energy eigenvalues of the 
strongly occupied orbitals reproduce the ionization po¬ 
tentials of small and large molecules [24, 25] accurately. 
We now turn our attention to the question of how the 
LROs compare with different sets of orbitals like the al¬ 
most exact NOs taken from MCSCF calculations, those 
from the full RDMFT minimization, and the KS orbitals 
of LDA. 


III. RESULTS AND DISCUSSION 

In this section, we perform full- and local-RDMFT cal¬ 
culations for the He atom, and the H 2 , molecule at 3 
different intermolecular distances, the equilibrium geom¬ 
etry with R = lA a.u., an intermediate distance with 
R = 3 a.u., and at the dissociation limit, R = 8 a.u., 
in order to examine the impact of the local potential 
constraint on the optimized orbitals. For full- and local- 
RDMFT we employ a few representative approximations, 
namely the Muller [2], BBC3 [5], Power [6, 14], ML [9], 
and Lowdin Schull (LSH) [27] functionals. The LSH func¬ 
tional is the exact one for the 2-electron singlet case in the 
full minimization. Apparently, in local-RDMFT, with 
the assumption of a local potential, LSH is no longer 
exact and not necessarily better than other approxima¬ 
tions. For comparison, we also calculated the NOs ob¬ 
tained within the complete active space MCSCF method. 
As these are very accurate approximations to the exact 
NOs, we will refer to them as “exact” NOs in the follow¬ 
ing. Full minimization orbitals from the LSH functional 
obviously coincide with those from the MSCF method, 
since we discuss only two electron singlet systems. Fi¬ 
nally, we compare the LROs with restricted HF (RHF) 
and LDA orbitals. All calculations were performed with 
the HIPPO computer code [29], except for MCSCF for 
which we used the GAMESS code [30]. In all calcula¬ 
tions, the cc-PVTZ Gaussian basis set was employed. For 
simplicity, and in order to avoid numerical noise intro¬ 
duced in weakly occupied orbitals, we compare orbitals 
with an occupation larger than 10“^. In order to com¬ 
pare orbitals qualitatively we plot them along the radial 
direction for the He atom and along the intermolecular 
axis for the H 2 molecule. To keep the figures uncluttered 
we plot only the orbital from one RDMFT functional 
(full and local), together with the MCSCF, the LDA and 


Distance 

BBC3 

LSH 

Exact 

(A) 

Local 

EKT 

Local 

EKT 


0.77 

15.98 

15.94 

15.98 

16.20 

16.20 

3.00 

12.75 

12.44 

12.89 

13.26 

13.26 


TABLE I. IPs (in eV) from the HOMO eigenvalues of local- 
BBC3 and local-LSH for H 2 at equilibrium and at 3 A 
compared with values from EKT with BBC3 and LSH full- 
RDMFT functionals. The exact values, which are obtained 
from the energy difference of the positive ion and the neutral 
system, are also included. 


the RHF ones. In addition, as a more quantitative com¬ 
parison, we calculate the overlaps of the corresponding 
orbitals obtained with full RDMFT and local-RDMFT, 
for all the functionals employed, with the exact NOs and 
KS-LDA orbitals. 

Before discussing the nature of the LROs, we present 
more evidence that static correlations can be correctly 
described by local-RDMFT if the corresponding IRDM 
functional has this property. In Fig. I, we show the oc¬ 
cupation numbers of the Itr and Icr* spin-orbitals of the 
H 2 molecule obtained with four different functionals in 
local-RDMFT. We see that the dependence of these oc¬ 
cupations on the distance is qualitatively correct going to 
the correct dissociation value of 1/2. Occupations from 
Local-LSH and local-BBC3, in particular, are quite accu¬ 
rate compared with the exact ones. The ML functional 
is not included in this plot as it does not yield a qualita¬ 
tively correct dissociation curve for H 2 in either the local 
or the full RDMFT version. As H 2 at large internuclear 
separation is a prototype correlated system, it is of par¬ 
ticular interest to see if the negative of the orbital energy 
of the HOMO is a good approximation for the ionization 
potential (IP). In the Table I, we show the calculated 
IPs at two internuclear separations, at equilibrium and 
at a distance of 3 A. Results from full RDMFT using 
the extended Koopman’s theorem (EKT) [26, 31-34] are 
also shown. EKT is proven to yield the IP exactly for 
the exact functional of the IRDM. It also provides a new 
orbital representation in terms of the canonical orbitals 
which can be associated with the electronic spectrum 
[34, 35]. Finally, in Table I, we also include the exact 
results within the particular basis set which are obtained 
using the definition of the IP i.e. the total energy differ¬ 
ence E{N — 1) — E{N), where E is the total energy and 
N is the number of electrons. As expected, the IP of the 
EKT method is exact for the LSH functional. Local-LSH 
and BBC3 functionals on the other hand are very accu¬ 
rate at the equilibrium distance and rather satisfactory 
at large separations as well, as shown in the Table 1. 

We continue with discussing the nature of the LROs. 
The helium Is orbital appears to be identical for all meth¬ 
ods considered in this work showing the typical shape of 
a Is orbital. This orbital is strongly occupied in local- 
and full-RDMET calculations and its occupation varies 
between 1.975 and 1.994. The He 2s orbital, on the other 
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hand, shown in Fig. 2, is strongly dependent on the cal¬ 
culation method. The occupation of this orbital varies 
between 0.0019 and 0.011. Within the local-RDMFT ap¬ 
proach, all four different functionals yield an almost iden¬ 
tical form on the scale of the plot, hence, only the one 
obtained with the Muller functional is included in Fig. 2. 
Also, the orbitals from the Muller and Power functionals 
in the full RDMFT calculation are identical, and again, 
only the one given by Muller is shown. Interestingly, all 
LROs from local-RDMFT are almost identical with the 
corresponding LDA orbital and are also close to the RHF 
one. Orbitals from full-RDMFT are closer to the exact 
NO with the orbital from BBC3 being the closest. 



local-RDMFT 

Exact NO LDA 

full-RDMFT 

Exact NO LDA 

Muller 

0.85803 

0.99964 

0.99657 

0.89302 

BBC3 

0.85748 

0.99967 

0.99956 

0.86825 

ML 

0.85878 

0.99953 

0.98462 

0.75565 

Power 

0.85792 

0.99963 

0.99620 

0.89507 

LSH 

0.85801 

0.99965 

0.99999 

0.85369 


TABLE II. Overlaps of the Helium 2s orbital employing dif¬ 
ferent functionals within local- or full-RDMFT with the cor¬ 
responding exact NO and LDA orbital. 


The overlaps of helium 2s orbital obtained with differ¬ 
ent methods are shown in Table II. The overlaps of the 
full-RDMFT orbitals with the exact NO are larger than 
0.98 for all approximations. On the contrary, the local- 
RDMFT orbitals for all functionals are almost identical 
to the LDA orbital with overlaps larger than 0.999. 

In Figs. 3(a) and 3(b) we plot the Icr and Icr* orbitals, 
for the hydrogen molecule at equilibrium, along the in- 
termolecular axis. The first orbital looks exactly the 
same for both local- and full-RDMFT and all approxi¬ 
mate functionals. It is also very close to the exact NO 
and RHF orbitals while the LDA gives a slightly differ¬ 
ent orbital which is, however, still close to the NOs. The 
overlaps between the RDMFT orbitals, both local and 
full, with the exact NO never deviate from one by more 
than 3 ■ 10“^. The overlaps with the LDA orbital are only 
marginally smaller with a deviation from one of about 
3 • 10“^. The antibonding Icr* orbital again looks exactly 
the same for all local-RDMFT functionals tested and al¬ 
most coincides with the corresponding LDA orbital. The 
RHF orbital is also close to them. The full RDMFT NOs 
are again close to the exact NO, with the BBC3 being the 
closest. The Power full-RDMFT orbital lies between the 
Muller and the BBC3 orbitals and is omitted from the 
plot for clarity. These findings are again confirmed by 
the overlaps between the different sets of orbitals which 
are given, for the Icr* orbital, in Table HI: The overlap 
between the local-RDMFT orbitals and the LDA orbital 
being around 0.999 in all cases and the full RDMFT or¬ 
bitals having an overlap with the exact NO larger than 
0.99. 



local-RDMFT 

Exact NO LDA 

full-RDMFT 

Exact NO LDA 

Muller 

0.74144 

0.99892 

0.99537 

0.76873 

BBC3 

0.73858 

0.99911 

0.99954 

0.71960 

ML 

0.73719 

0.99920 

0.99767 

0.75372 

Power 

0.74059 

0.99896 

0.99138 

0.64364 

LSH 

0.73759 

0.99917 

0.99999 

0.71255 


TABLE III. H 2 Icr* orbital at equilibrium: Overlaps between 
local- and full- RDMFT orbitals with exact NO and LDA 
orbitals for different RDMFT approximations. 



local-RDMFT 

Exact NO LDA 

full-RDMFT 

Exact NO LDA 

Muller 

0.80117 

0.99771 

0.93510 

0.88263 

BBC3 

0.79625 

0.99850 

0.99952 

0.79624 

ML 

0.79297 

0.99983 

0.95454 

0.89306 

Power 

0.79368 

0.99975 

0.98611 

0.75053 

LSH 

0.78874 

0.99991 

0.99999 

0.78008 


TABLE IV. H 2 2cr orbital at intermediate distance 3 a.u.: 
Overlaps between local- and full-RDMFT orbitals with exact 
NO and LDA orbitals for different RDMFT approximations. 

At an internuclear distance of 3 a.u. one finds practi¬ 
cally the same orbital for the Icr and Icr* orbitals for all 
methods considered in this work. Therefore, they are not 
discussed in more detail. In Fig. 4 we plot instead the 
2a orbital for this system, the only orbital, with a sig¬ 
nificant occupation, which differs depending on the em¬ 
ployed method. As in all other cases considered, local- 
RDMFT gives exactly the same orbital for all functionals. 
The corresponding LDA and RHF orbitals look similar 
and are close in shape to the local-RDMFT orbitals. On 



FIG. 4. H 2 2cr at internuclear distance 3 a.u.: All local- 
RDMFT functionals look almost the same. The full-RDMFT 
orbitals differ slightly from each other but are all closer to 
the exact natural orbital than to the LDA and RHF orbitals. 
Only the local and full-Miiller orbitals are shown. LDA and 
RHF orbitals almost coincide 
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1 ' 1 ' 1 

- Exact NO 



- LDA 



- - RHF 


- 

■ full-BBC3 

- 


■jA, ■ local-BBC3 A 

P \ 



_1_ 1 _1_ 1 _1_ 



_I_I_I_I_I_I_I_ 

-8.0 -4.0 0.0 4.0 8.0 


X (a.u.) 



FIG. 5. H 2 Icr (top) and H 2 Icr* (bottom) orbitals at disso¬ 
ciation. Essentially all functionals for full and local-RDMFT 
yield a symmetric and antisymmetric superposition of two 
hydrogenic orbitals except full-BBC3 that gives two localized 
orbitals instead as discussed in the text. 

the other hand, the orbitals obtained by full-RDMFT dif¬ 
fer significantly. While the BBC3 full RDMFT orbital is 
very similar to the exact NO, the Power, the ML and the 
Miiller orbitals differ from it especially close to the nu¬ 
clei. Overall, however, the full RDMFT orbitals are sig¬ 
nificantly closer to the exact NO than the local-RDMFT 
ones. This is also reflected in the overlaps given in Table 
IV which show a large overlap between the full RDMFT 
orbitals and the exact NO while the local-RDMFT or¬ 
bitals are rather closer to the LDA orbital than to the 
exact one. 

In Figs. 5 top and bottom panels, we show the first 
two orbitals for the hydrogen molecule at the dissocia¬ 
tion limit. Since most orbitals look almost identical on 
the scale of the plot we only included some examples to 
show the general behavior. At this limit the two orbitals 
plotted are the symmetric and antisymmetric combina¬ 
tion of hydrogen atom Is orbitals for all approximations 
except full-BBC3. The reason for this difference is be¬ 
cause this approximation is not invariant under a delo¬ 
calization unitary transformation in a subspace of local¬ 
ized NOs with degenerate occupation numbers. As it 
was shown in Ref. 36, BBC3 favors energetically orbital 
localization and, in that way, it remains size consistent. 
Local-BBC3, since it is based on the same energy func¬ 
tional, should also favor orbital localization. The reason 


Functional 

Alda 

AcLDA 

Muller 

-6.32 

-1.27 

BBC3 

-6.24 

-1.17 

ML 

-4.96 

-0.48 

Power 

-5.54 

-0.50 


TABLE V. The average differences (in eV) of the local- 
RDMFT eigenvalues, for the orbitals considered in the present 
work, from those of the LDA (Alda) and cLDA (Aclda). 

it does not is that the single particle Hamiltonian with 
the local potential is invariant in the subspaces of de¬ 
generacies (an occupation number degeneracy leads to a 
single particle energy degeneracy). Thus, for the effective 
Hamiltonian it is equivalent to get localized or delocal¬ 
ized pairs of orbitals and the choice is mostly a technical 
matter of the diagonalization routine. In our case, the 
result with local-BBC3 resembles that of the other ap¬ 
proximations, i.e. a pair of symmetric and antisymmet¬ 
ric delocalized orbitals. This small discrepancy between 
the full and the local version of a given functional can 
be resolved by a subsequent minimization of the original 
functional in the subspace of degeneracy. Generally, we 
should keep in mind that the local Hamiltonian is invari¬ 
ant under orbital rotations while the energy functional 
might not be. In such cases, a small scale additional 
minimization of the original functional in the degenerate 
subspace is required. 

Since for all systems studied here, the local-RDMFT 
orbitals are close to the LDA ones, one might wonder if 
the same is true for the corresponding energy eigenvalues. 
We found this not to be the case and the average differ¬ 
ences from the LDA eigenvalues for different functionals 
are of the order of 5-6 eV, for the orbitals considered 
in this work, as one can see in Table V. The relative 
differences for all functionals lie in the range of 32% to 
463%. This is largely expected since local-RDMFT en¬ 
ergy eigenvalues have been proven to give accurate ap¬ 
proximations to the experimental IPs [24, 25] contrary to 
the LDA ones. The big difference from the LDA eigen¬ 
values can be partly attributed to the constraints (3) 
and (4). If one compares the local-RDMFT eigenvalues 
with those from LDA calculations where the same con¬ 
straints are imposed [28] (referred to as constrained LDA 
(cLDA)) one finds an improved agreement in all cases, as 
shown in Table V. The average energy differences in this 
case are in the range of 0.5-1.3 eV for the different local- 
RDMFT approximations while the relative differences for 
all functionals are in the range of 0.2 - 64%. 

The close resemblance of local-RDMFT orbitals to 
LDA ones suggests that, as far as static correlations are 
concerned, one could completely avoid the orbital opti¬ 
mization. Instead one could perform a single minimiza¬ 
tion of occupation numbers using KS orbitals, for exam¬ 
ple from a LDA calculation, and still obtain a good de¬ 
scription of static correlation as local-RDMFT was found 
to provide. In order to obtain non-idempotent solutions 
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for the occupation numbers, it is essential to use 1-RDM 
approximations for the energy expressions and not den¬ 
sity ones [24]. The resulting method could provide an 
adequate description of strong static correlation at a com¬ 
putational cost which is only slightly larger than a density 
functional calculation since the only additional calcula¬ 
tion is a one-shot optimization of the occupation numbers 
which is rather cheap. In addition, a single evaluation of 
the effective potential of local-RDMFT with the modi¬ 
fied occupation numbers would offer an improved quasi¬ 
particle energy spectrum compared to the initial LDA 
calculation. 


IV. CONCLUSION 

In conclusion, this work focuses on the nature of the 
optimal orbitals, LROs, of local-RDMFT. This approach 
was introduced recently as a bridge between RDMFT and 
DFT. It is based on the idea of minimizing functionals of 
the 1-RDM under the additional condition that the LROs 
satisfy a single-particle Hamiltonian with a multiplicative 
KS-like potential. Apart from computational efficiency, 
the motivation for the development of local-RDMFT is 
that there exists a single-electron spectrum associated 
with LROs which was proven quite useful in describing 
spectral properties like ionization energies of molecular 
systems. 

According to the findings of this work, LROs from 
local-RDMFT are much closer to the Kohn-Sham or¬ 
bitals from LDA than the exact NOs. Orbitals from full- 
RDMFT on the other hand resemble the exact NOs much 
closer than LROs of local-RDMFT. This fact is a demon¬ 
stration of the hybrid nature of local-RDMFT: the com¬ 
bination of non-idempotency, i.e. fractional occupancies. 


with KS-like orbitals. Interestingly, the nice features of 
local-RDMFT, like the correct description of molecular 
dissociation and a single electron spectrum, are combined 
with KS-like orbitals due to the non-idempotency which 
is introduced through the adoption of functionals of the 
1-RDM. 

The present work can be seen as the initiative for the 
quest of approximations that combine density function¬ 
als leading to KS-like orbitals with approaches based on 
the 1-RDM that lead to fractional occupancies for these 
orbitals. Such schemes would combine the simplicity of 
DFT methods with advanced features of RDMFT, like 
the description of effects that require to go beyond the 
single determinant approximation. The similarity be¬ 
tween the local-RDMFT orbitals and KS-LDA orbitals 
suggests that combining a density functional calculation 
with a single optimization of the occupation numbers 
with a 1-RDM functional could provide an inexpensive 
way for the description of strongly correlated systems in 
a hybrid density/density matrix framework in the future. 
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